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ABSTRACT 

We use a generalised procedure for the combined likelihood analysis of different cosmological 
probes, the 'Hyper- Parameters' method, that allows freedom in the relative weights of the raw 
measurements. We perform a joint analysis of the cepheid-calibrated data from the Hubble 
Space Telescope Key Project and the baryon mass fraction in clusters to constrain the total 
matter density of the universe, £2 m , and the Hubble parameter, h. We compare the results 
obtained using Hyper-Parameters method with the estimates from standard \ 2 analysis. We 
assume that the universe is spatially flat, with a cosmological constant. We adopt the Big- 
Bang nucleosynthesis constraint for the baryon density, assuming the uncertainty is Gaussian 
distributed. Using this and the cluster baryon fraction data, we find that the matter density and 
the Hubble constant are correlated, Q. m h° 5 0.25, with preference for a very high h. To break 
the degeneracy, we add in the cepheid-calibrated data and find the best fit values (i2 m , h) = 
(0.26 



1-0.06 
-0.06' 



0.72jTo2) (68 per cent confindence limits) using the Hyper-Parameters approach. 



We use the derived Hyper-Parameters to 'grade' the 6 different data sets we analyse. Although 
our analysis is free of assumptions about the power spectrum of fluctuations, our results are in 
agreement with the A-Cold Dark Matter 'concordance' parameters derived from the Cosmic 
Microwave Background anisotropics combined with Supernovae la, redshift surveys and other 
probes. 

Key words: cosmology: observations - cosmology:theory - large-scale structure of universe 
- galaxies: clusters - methods statistical - X-ray: clusters. 



1 INTRODUCTION 

Combining different cosmological observations has become an es- 
sential and common approach in cosmology. A number of groups 
(e.g. Gawiser & Silk 1998; Webster et al. 1998; Lineweaver 1998; 
Eisenstein, Hu & Tegmark 1999; Efstathiou et al. 1999, 2002; Bri- 
dle et al. 1999, 2001a; Bahcall et al. 1999; Lahav et al. 2002) 
investigated a range of cosmological parameters by joint analysis 
of various cosmic probes, e.g. the Cosmic Microwave Background 
(CMB), Supernovae la (SNIa) and redshift surveys. 

It is well known that a simultaneous analysis of different probes 
is essential in finding tight constraints on the cosmological param- 
eters and breaking the intrinsic degeneracies inherent in any single 
measurement. While this is true, joint likelihood analyses pose sev- 
eral statistical problems. One of these problems arise when data 
sets are in disagreement (e.g. Press 1996). In this case, the general 
approach is to exclude the inconsistent measurements in a some- 
what ad-hoc way. A more objective approach to this problem was 
presented in Lahav et al. (2000) and Hobson et al. (2002), who 
generalised the conventional joint analysis by utilizing 'Hyper Pa- 
rameters' (hereafter HPs). The HPs provide a useful diagnostic in 
determining the relative weight that should be given to each exper- 



iment. Thus, this procedure gives an objective understanding as to 
which measurements are problematic and need further assessment 
of the random and systematic errors. The formalism for the HPs is 
given in Appendix B. 

In this paper we focus on Sl m , the ratio of matter density to 
the critical value necessary to close the Universe, and the Hubble 
parameter, Ho = WOh km s Mpc~'. We first consider the mea- 
surements separately and then combine them using the HPs method, 
as well as the conventional x 2 analysis. We assume a flat universe, 
n m + £1a = 1, in agreement with the latest CMB results (e.g. 
de Bernardis et al. 2002). We adopt the Big-Bang nucleosynthesis 
value for Clbh 2 , where Clb is the ratio of baryonic matter density to 
the critical density. 

We investigate the Hubble Space Telescope Key Project 
cepheid-calibrated data (Freedman et al. 2001) and estimates of 
gas mass fraction, / gas , in clusters of galaxies obtained from X-ray 
observations (Ettori & Fabian 1999; Mohr, Mathiesen & Evrard 
1999; Allen, Schmidt & Fabian 2002; following the earlier work of 
White et al. 1993). Unlike analyses (e.g. of the CMB and redshift 
surveys) that assume a power spectrum of fluctuations in a particu- 
lar scenario of dark matter, we have selected two probes that when 
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Sample 


h 


Error 


1.2 


36 Supernovae la (SNIa) 


0.71 


± 0.02 r ± 0.06 s 




21 Tully-Fisher Clusters (TF) 


0.71 


±0.03 r ± 0.07 s 


1.0 


1 1 Fundamental Plane Clusters (FP) 


0.82 


it 0.06 r ± 0.09 s 




Surface Brightness Fluctuations (SBF) for 6 Clusters 


0.70 


± 0.05 r ±0.06 s 


> 0.8 



SBF, SNIa, TF, FP: 



SBF, SNIa, TF, FP: Joint x' and HPs 



Table 1. The values of h and the 1-sigma random(r) and systematic(s) 
uncertainties for Cepheid-calibrated samples (from F01, and the references 
therein) 



combined measure f2 m with minimal assumptions, independent of 
the nature of dark matter. The main assumptions we made in our 
analysis are as follows: (i) the local Ho measurements are typical of 
the entire universe, (ii) the clusters of galaxies are representitive of 
the matter distribution on large scales (but see e.g. Bahcall & Com- 
erford 2002), (iii) all the systematic errors of the measurements are 
included. 

In Sections 2 and 3, we introduce the data and discuss the 
methods used to produce likelihoods for the individual data sets. 
We combine the two data sets in Section 4. In Section 5, we discuss 
our results and compare them with other studies of baryon fraction 
and independent measurements from the CMB and the 2dF galaxy 
redshift survey. 



2 HUBBLE CONSTANT FROM CEPHEID-CALIBRATED 
DISTANCES 

One of the most important results on the Hubble constant comes 
from the Hubble Space Telescope Key Project (Freedman et ah 
2001, hereafter F01). The group has used the cepheid period- 
luminosity relations to obtain distances to 31 galaxies and cali- 
brated a number of secondary distance indicators measured over 
distances of 400 to 600 Mpc. The values of h derived using these 
methods and the uncertainties (random (r) and systematic (s)) are 
summarised in Table 1. Combining the measurements in Table 1 by 
several statistical methods they derive the final result as h = 0.72± 
0.03 r ± 0.07 s . Given the importance of this widely quoted result, 
we perform a more principled statistical analysis than used by F01 
to test robustness of this value. 

We use the raw data given in the tables in F01 for Surface 
Brightness Fluctuations (SBF), SNIa, Tully-Fisher (TF) and Fun- 
damental Plane (FP). We combine these using the standard joint x 2 
and then the HPs approach. Each method is affected by both system- 
atic errors which are common to all of the methods (e.g. the adopted 
distance modulus to the Large Magellanic Cloud, metallicity depen- 
dence of cepheid period-luminosity relation and reddening by dust) 
and systematic errors which are specific to each method. In princi- 
ple, some of systematic errors can be modelled and incorporated in 
the x analysis by adding extra parameters and marginalising over 
them (for an application to CMB data see Bridle et al. 2002). How- 
ever, we follow F01 and for the variance, a 2 , in the x 2 analysis we 
use the quoted random and systematic errors (from Table 1), added 
in quadrature. We also test the results when we assume random 
errors only (see Table 3). In both cases we find the best fit h = 0.72 
with the HPs method. This value is in very good agreement with 
the HST Key project result and also with the standard x~ analysis 
we performed (Table 2). 

The HP values we obtain for each case are given in Table 3. 
As can be seen from the HPs formalism in Appendix B, the HPs 
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Figure 1. Probability functions for the Hubble constant. The left plot shows 
the x 2 statistic for four Cepheid-calibrated distance indicators from F01: 
SBF(dashed-dotted line), SNIa(long-dashed),TF(dotted) and FP(dashed). 
The right plot shows the probabilities based on the joint x 2 (long-dashed) 
and the HPs approach (dotted) for the Hubble constant given the same data, 
with peak values of h = 0.72 and h = 0.73, respectively. See also Lahav 
2001b. 



can be interpreted as either an indication of misfit of the data set 
and the model (e.g. due to systematic problems with the data or 
to the wrong model), or as a rescaling of the quoted error bars 
(to a/a l/2 ). The HPs scores rank the quality of the methods as 
follows: SBF (the highest), TF, SNIa and FP. This is in accord with 
comments on the precision and systematics of each method given 
by F01. The FP measurements alone give the most discrepant value 
of h (Figure [l] and Table 2), and they got the lowest HP (even 
when the systematic errors were not included in the analysis; see 
Table 3). This is an example where HPs flag a problem with the 
measurements. Although the discrepancy of the h value obtained 
using FP data alone does not affect the final value derived by joint 
analyses of FP, SNIa, TF and SBF, one needs to be careful about 
using the standard x~ ■ A X 2 analysis of discrepant results may lead 
to a wrong joint answer. A clear example of this problem arises if 
we use only the FP and the TF data, which have the most similar 
scatter (see Figure |l|) and x 2 values (see Table 2) out of all the 
probes. We find that the joint likelihood function for h peaks at 0.76 
using HPs analysis, whereas the best value of h is 0.80 using the x 2 
analysis. In other words, HPs analysis grades the TF measurements 
as more reliable than the FP data and hence gives the joint value of 
h closer to the best value obtained using TF measurements. Using 
standard joint analysis, the joint likelihood function peaks between 
the likelihood functions of the two probes, as expected. It is also 
interesting to note that indeed F01 attached to the FP data the largest 
systematic errors (see Table 1). 

As can be seen in Figure [j], the joint probability functions 
obtained using the two approaches (x 2 and HPs) are both dominated 
by SNIa data, since SNIa as an indicator possesses the smallest 
random scatter. 



3 X-RAY OBSERVATIONS OF CLUSTERS OF GALAXIES 

Clusters of galaxies provide various means for the study of the cos- 
mological parameters. Clusters are the largest gravitationally bound 
structures in the Universe, and therefore generally are assumed to be 
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Sample 


9 

X' 


Best h 68 


per cent confidence limits 


36 SNIa 


19.0 


0.72 


0.70 < h < 0.73 




7.6 


0.74 


0.70 < h < 0.76 


11 FP 


4.9 


0.88 


0.84 < h < 0.92 


6SBF 


1.7 


0.72 


0.66 < h < 0.76 


Joint x 2 (f° r 74 data points) 


57.0 


0.73 


0.70 < h < 0.77 


Table 2. The conventional x 2 


analysis using the raw F01 data. For each 


sample the best fit value of h and the 


X 2 value at this point are given. 


Sample 




HPs 


HPs 






(r only) 


(r+s) 


36 SNIa 




0.3 


1.9 


21 TF 




1.9 


2.7 


11 FP 




0.2 


0.5 


6SBF 




2.8 


3.4 


Best h 




0.72 


0.72 


68 per cent confidence limits 


0.69 


< h < 0.76 


0.68 < h < 0.77 



Table 3. The Hyper-parameters analysis using the raw F01 data. For each 
case the HP is given. The first column is shows the HPs obtained using 
random errors only, while the second column is for random+systematic 
errors, added in quadrature. 



the tracers of matter distribution on large scales. Provided that £ib/T, 
the baryon density, can be inferred from primordial nucleosynthesis 
abundance of the light elements, the cluster baryon fraction, fb = 
can then be used to constrain £l m and h (White et al. 1993, 
Steigman, Hata & Felten 1999, Ettori 2001). The baryons in clus- 
ters are primarily in the form of X-ray emitting gas that falls into the 
cluster halo and secondarily in the form of stellar baryonic mass. 
Hence the baryon fraction in clusters is estimated to be 



fb — ~ /gas + /gal, 



(1) 

M g ai/M gr av and 



where f b = M b /M gmv , / gas = M gx /M grM , / gaJ 
M glav is the total gravitating mass. 

We consider two different datasets of gas mass fraction esti- 
mates. The first was published in Ettori and Fabian ( 1 99W (hereafter 
EF99). This is a sample of 36 relaxed galaxy clustersEI with high 
X-ray luminosities (Lx il 10 45 erg s _1 ) and a redshift range from 
0.05 to 0.44. The second data set is taken from Mohr, Mathiesen 
and Evrard ( 1 999) (hereafter MME99) . This is an X-ray flux-limited 
sample of 45 clusters with redshifts between 0.01 and 0.18. Both 
groups use ROSAT PSPC surface brightness profiles and interclus- 
ter medium (ICM) temperatures from ASCA, Ginga and Einstein 
MPC observations for their analyses. The gas mass fraction esti- 
mates, /g as , in both datasets are derived assuming the gas is isother- 
mal and in hydrostatic equilibrium within the limiting radius, R A , 
where A is the mean overdensity of the total mass in a cluster relative 
to the background value. 



In our analysis, we use 35 clusters instead of 36. Cluster A3888 is ex- 
cluded from the sample since this cluster's gas temperature was not obtained 
from X-ray observations but from the gas temperature-optical velocity dis- 
persion relation. 



To allow for the variation from cluster to cluster, for the present 
analysis we use the individual / gas per cluster, rather than a global 
average. As the clusters are at relatively high redshifts, / gas is a 
function of £2 m and h (see Figure 2).This is due to two cosmological 
effects (see discussion in, e.g. Ettori 2001). 

The first dependence is due to the angular diameter distance, 
dang, which is a function of h and £2 m in a flat universe with cosmo- 
logical constant. For an isothermal gas in hydrostatic equilibrium, 
/ g as can be calculated through the surface brightness profile and gas 
temperature (cf. Ettori & Fabian 1999; Mohr, Mathiesen & Evrard 
1999). If the surface brightness profile is measured as the integra- 
tion of the thermal bremsstrahlung emissivity of the ICM along the 
line of sight, one finds that / ga8 is proportional to dl n l (see A ppen dix 
A). We calculate the angular diameter distances using Eq. <\4 and 
vary h and f2 m . We find that / gas decreases significantly (about 40 
per cent) in a low density universe with h > 0.5. 

The second cosmological dependence is weaker. The hydro- 
dynamical simulations (see e.g. Evrard et al. 1996) show that for an 
Einstein-de Sitter cosmology, a mean overdensity, A = 500 defines 
a region within where the assumptions of isothermal gas in hy- 
drostatic equilibrium are valid. Therefore, in both data sets, / gas has 
been estimated at A = 500. However, A is function of the cosmo- 
logical parameters and increases as Q. m decreases in a universe with 
positive £^a- -Ra is proportional to (f2 m A) -0 ' 5 (Eq. kl in Appendix 



A) and combining this with the radial dependence of the gas mass 
fraction near R A (/ gas oc r ' 2 ), we get / gas oc (fl m A) -01 . Figure ^ 
illustrates the dependence of / gas on h and £2 m . We see that / gas is 
mainly sensitive to h. 

We also include the stellar contribution in cluster galaxies, 
/gai, to the baryon fraction estimate. We use a global correction 
/gai = (0.01 ± 0.005)/i _1 (Fukugita et al. 1998), which is based 
on mass-to-light ratio, M/L — 4.5 ± 1 per galaxy, after averaging 
different galaxy types, with no dependence on h, hence the h- 
dependence in / ga i comes solely from the h dependence of Mg rav 
on the cluster distances. We note that some estimates of M/L 
per galaxy depend on the distances to galaxies, and hence on h, 
resulting in /i-independent / ga i (e.g. Wang et al. 1999 quote / ga i = 
0.013). Ideally, / ga i should be estimated per cluster, but these data are 
currently unavailable, and in any case it is only a small contribution 
to the total baryon fraction. 

We vary the matter density of the universe in units of the crit- 
ical density, Q. m , and the present day value of the reduced Hubble 
constant, h, over the parameter spaces [0.01, 0.5] and [0.5, 2.5], re- 
spectively. We also put a prior on fib h 2 so that £lbh 2 = 0.019±0.002, 
95% confidence levels, (e.g. Buries et al. 2001) and marginalise over 
the uncertainty range in order to get the likelihoods of f2 m and h. 
We note that fixing Clbh 2 to 0.019 tightens the constrains in the f2 m , 
h plane but does not change the best fit points significantly. 

We compute the x 2 s using Eq. B2 in Appendix B: 



(2) 



where the sum is over the number of clusters, the baryon fraction is 
given by / b = / gas + /gai and w b = Clbh 2 . 

To obtain a qualititive understanding we consider Eq. [|and use 
the first cosmological dependence that dominates over the others to 
an approximate relation between fl m and h: 



On 



u>b h 



LOb h 2 

/gas + /gal ^ 0.0Sh- t S +0.01/1- 



(3) 



© 0000 RAS, MNRAS 000, 000-000 



4 Pirin Erdogdu, Stefano Ettori, Ofer Lahav 



EF99 MME99 




Gas Fraction f g „ Gas Fraction f fi m !! m 



Figure 2. The histograms of the gas mass fraction /b observed in 35 clusters 
(EF99). The left panel shows the dependence on h for fixed Q. m = 0.3: 
h = 0.5 (solid line), h = 0.75 (dashed line) and h = 1.0 (thick solid line). 
The right panel shows the dependence on £2 ra for fixed h = 0.75: Q. m = 0.1 
(solid line), £l m = 0.5 (dashed line) and £2 m = 1.0 (thick solid line). 



For the observed distribution of the gas mass fraction of about 0.08 
h~ 1-5 (e.g. EF99) and reasonable values of Hubble constant, the 
above relation can be estimated as Sl m h 05 w 0.25. 

The two datasets do have some clusters in common, however, 
we see that our conclusions are not affected if we consider the sam- 
ples separately. The results are plotted in Figure ^. It is seen that 
in both data sets the value of h is unacceptably high at the 68% 
confidence level (but still lower than Hubble's original value of 500 
km/sec/Mpc!) and the best fit values of Sl m are relatively low. The 
parameter values (Sim, h) at the best fit points (with 68 per cent con- 
fidence limits) are (0.11±^,l-73±o;||) and (0.18iJ^,l.ir!:g;}g) 
for EF99 and MME99 data, respectively. The standard joint analy- 
sis of the data sets yields (Sim, h) = (0.171^,1.23^*) and with 
HPs method we obtain (Sim, h) = (0.15±° °£, 1.33±£^). It is easy 
to see from Figure ^| that both separate and joint analyses of the 
data imply that Sl m = 1 is ruled out at very high confidence level. 
The standard x~ and HP analyses, shown on bottom left and right 
panels in Figure | yield slightly different results. The HP values are 
0.6 for the EF99 and 0.1 for the MME99 sample. The low weight 
given to the MME99 data may indicate possible systematic errors, 
an underestimation of the random errors or incomplete modelling. 
The HPs obtained suggest that the EF99 sample is more reliable 
than the MME99 data. The galaxy clusters in the EF99 sample were 
selected for their high X-ray luminosity and relaxed morphology, 
whereas the selection criteria in MME99 was to build a flux-limited 
sample. Therefore, we see at least two effects that can make the 
EF99 sample more robust in providing a stable central value of 
the gas fraction: (i) the systematics due to non-homogeneous ob- 
jects are more under control and (ii) the observed dependence of 
the value of gas mass fraction upon the plasma temperature (and 
luminosity, as consequence of the tight L — T relation observed in 
galaxy clusters; e.g. Ettori, Allen & Fabian 2001 for an application 
of HPs to this relation) makes the selection in luminosity (instead 
of flux) a more robust way to select objects in the upper end of the 
gas mass fraction distribution. 

Being concerned about the high Hubble constant, we applied 
the Bootstrap method (Efron 1982) to baryon mass fraction data to 
ensure that there are no outlying clusters which could alter the sig- 



Figure 3. The X-ray Cluster Likelihood Functions in the {£2 m ,/i}-plane, 
after marginalisation over Cl^h 2 . The results are compared to estimates 
from the 2dF Galaxy Redshfit Survey: Q. m h ss 0.2 (dashed lines; Percival et 
al. 2001) and from the CMB: Q m h 2 ps 0.15 (dotted lines; Netterfield et al. 
2002 and Pryke et al. 2002) and the age of the universe to as 14 Gyr (solid 
lines; Knox et al. 2001). The contours denote the 68, 95 and 99 percent 
confidence regions. 

nificance of the obtained best fit values. We created 2000 synthetic 
catalogs selected from the two samples. The histograms of the best 
fit points for h and Sl m for these bootstrap realizations are in very 
good agreement with the presented results. 

The joint analysis of cosmological probes (e.g. Bahcall et al. 
1999 and Efstathiou et al. 1999) suggests a flat Universe with 
Sl m ~ 0.3. We fix Sim to 0.3 and to 0.2 and plot 1 dimensional 
likelihood distributions of h. These plots are shown in Figure ^. 
For fixed £2 m = 0.3 the best fit points of h are 0.64 and 0.78 for 
EF99 and MME99, respectively. The joint 1 -dimensional likelihood 
distribution peaks at h = 0.72 using the standard and at h = 0.66 
using the HPs approach. Note the good agreement of the best fit 
value of h for both data with the F01 result. When £2 m is fixed to 
be 0.2, the plots shift significantly to the right, with best fit points 
h = 1.03, 1.11, 1.07 and 1.03 for EF99, MME99, standard joint 
analysis and HPs approach, respectively. 

We conclude that the baryon fraction data on its own cannot 
constrain each of the two parameters, but only their combination, 
Slmh ~ 0.25. Therefore, although the likelihood peak is at a 
high value of h and a low value of Sl m , we should not attach much 
significance to the individual values. To constrain the individual 
values of h and £2 m we now combine the baryon fraction data with 
the cepheid sample. 



4 COMBINING THE BARYON FRACTION AND THE 
CEPHEID DATA 

We present our results for the combined analysis in Figure |5| It can 
be seen that the high contour regions obtained for h with the cluster 
data have decreased significantly. The confidence regions are a lot 
tighter than the confidence regions of the single data sets alone, thus 
giving stronger constraints. 

Table 3 and Table 4 show the best fit values and 68% confi- 
dence limits of h and £2 m , respectively. The best fit points of the 
parameters lie within the confidence limits, indicating that the like- 
lihood distributions are well behaved. The dominating data sets are 
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Figure 4. The 1-d likelihood distributions of the Hubble Constant for fixed 
Q. m = 0.3 (dashed lines) and £2 m = 0.2 (solid lines) and marginalised over 
Q b h 2 . Theft gure illustrates that h is highly sensitive to the assumed H 

m, as 

h and £2 m are highly correlated. 



Figure 5. The likelihood functions from combining the Cluster data with 
Cepheids.The results are compared to estimates from the 2dF Galaxy Red- 
shfit Survey: Q. m h sa 0.2 (dashed lines) and from the CMB: Cl m h 2 PS 0.15 
(dotted lines) and the age of the universe to ~ 14 Gyr (solid lines), as in 
Figure 3. The contours denote the 68, 95 and 99 percent confidence regions. 



Analysis Best fit value 68 per cent confidence limits 

Joint x 2 0.74 0.70 < h < 0.78 

HPs 0.72 0.70 < h < 0.76 



Table 4. The values derived for h using different analysis techniques. The 
68% confidence limits are given, calculated for each analysis by marginal- 
ising the likelihood function over Q, b h 2 and Cl m . 



the cepheid-calibrated distance indicators. This is expected, since 
the uncertainties on h are much larger for the cluster samples. 

The resulting HPs of our analysis are 0.6 (EF99), 0.1 
(MME99), 1.8 (SNIa), 2.7 (TF), 0.5 (FP) and 3.3 (SBF). The HPs 
are actually almost identical to the HPs derived for the individual 
data sets in section 2 and 3. We also see that the dominance of the 
cepheid data relative to the baryon fraction data. Since the cepheid 
data only constrain h, this leads to a narrower error bar for h in the 
HPs analysis compared with the joint x (which give equal weight 
to each of the 6 data sets). Since the cepheid data have no informa- 
tion about S2 m but have higher HP values, the error bar on £2 m is 
wider in the HPs analysis. 

Another interesting result is that the HPs analysis probes a 
parameter space in good agreement with the results from 2dF Galaxy 
Survey and CMB anisotropy measurements. The right panel of 
Figure ^ shows that our combined analysis is consistent with other 
measurements. 



Analysis 


Best fit value 


68 per cent confidence limits 


Joint x 2 


0.25 


0.22 < Q m < 0.29 


HPs 


0.26 


0.20 < On, < 0.32 



Table 5. The values derived for I2 m using different analysis techniques. The 
68% confidence limits are given, calculated for each analysis by marginal- 
ising the likelihood function over Q.\,h 2 and Cl m . 



5 DISCUSSION 

We have presented applications of a generalised procedure, 'Hyper- 
Parameters' (HPs), for analysing a set of different measurements. 
We performed a combined analysis of baryon fraction in clusters and 
cepheid-calibrated distances. We used the HPs formalism for joint 
analyses of the data sets to constrain the cosmological parameters 
h and £2 m (assuming a flat universe) and to check the reliability of 
the data sets. 

Using the baryon mass fraction in clusters, we obtained a very 
high best fit for h . However, Figure 3 shows a strong correlation 
between h with Cl m , Q, m h 5 ~ 0.25. The addition of the cepheid 
sample to the cluster data changed significantly the best fit values, 
h — 0.72 and CI = 0.26. This is not surprising, as the accurate 
cepheid-calibrated distances dominate the joint likelihood. 

Recently, Douspis et. al (2001) combined CMB data with a 
fixed value for baryon fraction fi, — 0.048/i~ 1,5 + 0.014 and found 
Cl m ~ 0.4, h ~ 0.6. Their results are in agreement with our derived 
combination, £l m h 5 « 0.25 and in marginal agreement with our 
best fit values for £l m and h, derived from the joint analysis of 
baryon fraction and the cepheid data. Our analysis differs from 
theirs in the sense that we have used individual gas fraction values 
for each of the 80 clusters in our sample, allowing self-consistency 
for the dependence of the gas fraction on Cl m . 

More recently, Allen et. al (2002) used gas mass fraction mea- 
surements of 6 relaxed clusters observed with the Chandra Obser- 
vatory, the HST Key Project value for Ho and the and Big-Bang 
nucleosynthesis value for fib h~ to constrain £2 m and fl A . The value 
they obtain for £2 m is again in very good agreement with our results. 

In our analysis, we have found that HPs assigned to the cluster 
samples are lower than the values assigned to cepheid-calibrated 
distances. Indeed, it is worth noticing that some systematic uncer- 
tainties affect / gas when it is estimated under the assumption of 
an isothermal gas in hydrostatic equilibrium with the dark mat- 
ter gravitational potential, like in this paper. Some aspects of the 
physics of the ICM still need further understanding (e.g. the mea- 
surement of the total mass changes in presence of a gradient in 
the temperature profile and/or supporting pressure from non ther- 
mal component) and will be addressed in the near future with the 
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current X-ray missions Chandra (Weisskopf etal. 2000) and XMM- 
Newton (Hasinger et al. 2001). Moreover, we have not taken into 
account all contribution from baryons in dark matter to the cluster 
baryon budget. 

There are also significant systematic uncertainties in the HST 
Key Project Result which should be addressed in the future. These 
uncertainties, discussed in detail in F01, are mainly due to the 
errors in the adopted distance modulus to the Large Magellanic 
Cloud (LMC) upon which the Cepheid period-luminosity rela- 
tion is strongly dependent. Furthermore, the effects of reddening 
and metallicity on the Cepheid period-luminosity relation maybe 
stronger than expected (Shanks et al. 2002). It is also worth noting 
that the local variations in the expansion rate due to large scale 
velocities may effect the accurate determination of Ho (eg. Turner, 
Cen & Ostriker 1992). 

Although our analysis is free of assumptions about the power 
spectrum of fluctuations, the results we obtain are in remarkable 
agreement with the A-Cold Dark Matter 'concordance' parameters 
(e.g Figure 5) derived from the Cosmic Microwave Background 
anisotropies combined with Supernovae la, 2dF galaxy redshift 
survey and other probes. 
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APPENDIX A: ESTIMATION OF THE GAS FRACTION 

For an isothermal plasma in hydrostatic equilibrium, the gas mass 
fraction, / gas , is calculated as the ratio of the mass of the gas and 
of the total gravitating mass, M grav , within the radius, Ra, at which 
the given mean overdensity of the total mass within cluster, A, with 
respect to the background value, £2 ra p c , is reached: 



M m (<R A ) _ 4tt J q R& p gas (r)r 2 dr 



/gas — 



M grav (< -Ra) kT^Ri 



Gfin 



( d\np m \ 
[ e In r ) 



(Al) 



where R A = 9 d allg is the physical radius, 9 is the angular separation, 
dang is the angular diameter distance, T gas is the ICM temperature, k 
is the Boltzmann constant, /i is the mean molecular weight in a.m.u. 
(~ 0.6), G is the gravitational constant and m p is the proton mass. 

The surface brightness, S(9) is given by the integral along the 
line of sight: 



s(9) « I /-:;,„/::;.//. 



(A2) 



Hence, the gas density, /9 gas is pro port ional to d ang ' 5 . Combining this 
with the other dependence in Eq 

"■■arm 



'porlii 

m 



/gas OC 



d anB oc h 



-<<ang 

For £2k = 0, d ana can be written as 



dang — 



H (l + z) ' I Jo |- 



(A3) 



if £2 m = 
otherwise (™) 
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The mean overdensity, A, for £\= is given by (e.g. Kitayama 
& Suto 1996; Henry 2000) 

A(O m , z) = A(£2 m = 1) ■ 1 + 0.4093 (fl m ; z - l) 9052 , (A5) 

where (n~' z - 1) = (£tf - 1)/(1 + z)\ 
The observed mass profile is calculated as 

M grav (< Ra) = (4/3)^1^.(1 + z) 3 A (A6) 

and is proportional to R& for a virialised system so that 

R A <x (fi m A)- ' 5 . (A7) 



-2 In P(w\D A ,D B ) =NaHxa) + NbHx\ 



(B9) 



To find the best fit parameters w requires us to minimise the above 
probability in the w space. It is as easy to calculate this statistic as 
the standard ^ 2 , and it can be generalized for any number of data 
sets. 

Since a and /3 have been eliminated from the analysis by 
marginalisation they do not have particular values that can be 
quoted. Rather, each value of a and (5 has been considered and 
weighted according to the probability of the data given the model. 



It can be shown that the 'weights' are a e ff 
both evaluated at the joint peak. 



and /3 c ff 



APPENDIX B: HYPER-PARAMETERS 

Lahav et al. (2000) (see also Bridle 2000, Lahav 2001a & 2001b), 
generalised the standard procedure of combining likelihoods using 
the 'Hyper-Parameters' approach. The conventional way of com- 
bining likelihood functions of different data sets is either to give all 
the sets the same statistical weight or to assign weights in an ad-hoc 
way. The 'Hyper-Parameters' method generalises this approach by 
assigning each data a relative weight. 

Given two independent data sets Da and Db (with Na and 
Nb data points respectively), our approach is to combine the x 2 s 
in the following manner: 



xfotoi = a X~A + P Xb , 



(Bl) 



where a and (3 are 'Hyper-Parameters' (HPs). The maximum likeli- 
hood of a given model is estimated by minimising the above quan- 
tity. We calculate the \ s using the equation below: 



X 



(B2) 



where the sum is over the number of measurements, Oi is the error 
for each data point and w is the vector of free parameters we wish 
to determine (e.g. £2 m and h). 

The HPs are eliminated by marginalisation over a and f3: 



P{w\D A ,D B ) = 



P(w,a,P\D A ,D B ) dad/3 



(B3) 



In order to evaluate the above integral, we use Bayes' theorem 
to write the following relations: 

P(w,a,f>\D A ,D B ) = nD^Ds^P) P(w,a,/3) ^ ^ 



P(D A ,D B ) 



and 



(B5) 



P(w,a,0) = P(w\a,0) P(a,P). 
We also assume the following: 

P(D A ,D B \w,a,p) =P(D A \w,a) P(D B \w,/3) , (B6) 



P(w\a,f3) — const. , 
and 

P(a,f3)=P(a) P(p) 



(B7) 



(B8) 



We take the prior probabilities in Eq. |B8| as Jeffreys' uniform priors 
in the log, P(lna) = P(ln/3) = 1. Assuming Gaussianity, we 
write P(-Da|w,o) oc ol Na ^ exp( — §Xa) an d similarly for 
Db - It then follows that the probability for the parameters w given 
the data sets is: 
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